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The non-resonant electromagnetic transmission of a normal-incident plane wave through a single 
hole in a perfect conductor metal slab of finite width is studied. The cases of rectangular and 
circular holes are treated in detail. For holes in the extreme subwavelength regime, in a film of finite 
thickness, the transmittance is shown to have the Rayleigh dependency upon the wavelength and, 
in addition, is mainly suppressed due to attenuation of the fundamental waveguide mode. In the 
limit of an infinitesimally thin screen Bethe's result is recovered for the circular hole. The numerical 
computations are fitted, providing expressions for the transmission in a wide region of parameters. 
We reformulate our results in terms of multipole expansion, interpreting the waveguide modes inside 
the hole as induced multipole moments. This result provides the link between the modal expansion 
method and the one based on a multipole expansion. 
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I. INTRODUCTION 

Electromagnetic (EM) wave transmission through 
apertures in perfect metal screens has been the subject 
of multiple studies. For a long time, up to the middle 
of XX century, theoretical treatments of diffraction by 
either opaque or metal bodies based on the Kirchhoff 
approach. This method consists in setting the fields 
on the body equal to their incident values. Bethe was 
the first to consider the diffraction by a small circular 
aperture of radius a in an infinitesimally thin perfectly 
conducting screen, providing rigorous analysis of the 
Maxwell's equations with the exact boundary condi- 
tions in 1944.12 He found analytically the pre-factor 
appearing in the well-known Rayleigh scattering de- 
pendency ~ (a/A)^ for the scattering cross-section of 
small objects in the long- wave length limit (a ^ A). 
Later on, Bowkamp improved Bethe's result, provid- 
ing additional terms in a series expansion of the trans- 
mission over a/ A. Since then, this kind of diffraction 
problem has become a classical chapter in many mono- 
graphs on clectromagnetism, e.g., Refs. [2131 

The interest in EM transmission through apertures 
has been renewed thanks to the discovery in 1998 of 
the enhanced optical transmission through an array 
of small holes.** A great deal of research has been de- 
voted to the transmission through periodical arrays of 
the holes. Nevertheless, the analysis of the diffraction 
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by a single hole in a film of finite thickness is still in- 
complete, as results found in the literature are only for 
fixed geometrical parameters. For example, Refs. I5I6I 
provide computed transmission spectra through a cir- 
cular hole in a film of a finite thickness in a perfect elec- 
tric conductor (PEC) slab. The case of a single circular 
hole in a real metal has been considered in Refs. I7I8I9I 
However, it is problematic to extrapolate these results 
to other parameters and other hole shapes. An at- 
tempt to represent the solution in an analytical form 
was undertaken in Ref. 1101 These authors derived the 
normalized cross section for a circular aperture (of ra- 
dius up to half of the wavelength in PEC) assuming 
that the magnetic current is uniform within the aper- 
ture. Nevertheless, this model shows poor quantitative 
agreement with both the Bethe-Bouwkamp results and 
the strict numerical calculations of Roberts .'^ More- 
over, it is limited to a screen with zero thickness. 

In this paper we study the optical transmission 
through single holes in PEC. We present new ana- 
lytical results, valid for a wide range of geometrical 
parameters, and provide a link between the modal ex- 
pansion method and the one based on multipoles. The 
paper is organized as follows: In Section II we de- 
scribe the modal expansion technique used for studying 
the transmission through the hole in a perfect electric 
conductor film of arbitrary thickness. In Section HI 
we provide analytical expressions for the transmission 
through holes of both rectangular and circular shapes 
in extreme subwavelength limit. We test our approach 
by applying it to the most unfavorable conditions for 
the method used (when the film thickness is zero) , and 
obtain an excellent agrement with known results. In 
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Section IV the square hole of moderate size is treated, 
when the wavelength is still larger than the hol e cutoff 
(therefore, the resonances found close to cutofP^ISIill 
will not be discussed there, as have already been ad- 
dressed befor^^^E^. Finally, in Section V, we make the 
link between the induced multipoles and the waveguide 
modes inside the hole. We discuss the importance of 
the fundamental waveguide mode both for the forma- 
tion of the induced dipole moments on both faces of 
the hole and for the coupling between these moments. 
In Appendix we explain the simplifications of Green's 
tensor for the small hole limit. 



II. THEORETICAL BACKGROUND 
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FIG. 1: (Color online) Geometry of the problem. 

In this section we briefly outli ne t he modal expan- 
sion formalism for the EM field ."^^EH Consider an EM 
wave incident onto a PEC film of finite thickness h con- 
taining a hole, see Fig{l] We first assume that both the 
bounding media and the medium inside the aperture 
is vacuum, but later on we generalize the results to ar- 
bitrary values of the optical indices of the semi-infinite 
media. 

For the representation of the EM field in-plane com- 
ponents we use Dirac's notations as follows. The pro- 
jection of a given Dirac's vector onto the position vec- 
tor |rt), with Tt — {x,y), simply yields the value of 
Dirac's vector in the position rt coordinate dependent 
field in the coordinate space, for example: 



(1) 



The projection of a vector onto another one is given 
by the scalar product 



{a\K) 



drt{a\rt){rt\Ky 



(2) 



means complex conjunction. 
Let us write in Dirac's notations the tangential com- 
ponents of the fields in the lower and upper half-spaces, 
expanding them over the continuum of plane waves 
(see the geometry in Fig. [T]): 



(z-h). 



(3) 



Here k = (k, a) represents both the in-plane wavevec- 
tor component k = (kx, ky) and the polarization, a = p 
or CT = s. The coordinate representation of the modes 
in the vacuum half-spaces reads 



(rt|k,s> 



ikr. 



|k,p) 



(4) 



The plane wave propagation constant is k^ — 
\/ — k^ with g = 27r/A. The summation operator in 
Eq. ^ includes both the integration over the k con- 
tinuum spectrum and the summation over the polar- 
izations: J2k. ~ (1/2'"')^ / '^^^ Inside the hole, we 
expand the tangential components of the field over the 
modes \a) of the infinite waveguide 

|E(z)> - E (^ae^9-^ + B„e-'«-^) (5) 

a 

where and are the amplitudes of the waveguide 
modes propagating (or decaying) forwardly and back- 
wardly with respect to z-axis direction; q^a represents 
the propagation constant of a waveguide mode with 
the label a. This label contains both the polarization 
of the waveguide mode and a "spatial" index related 
to the number of nodes of the field inside the hole. 

By matching the EM fields at the interfaces, and 
using the orthogonality of the modes, we arrive at a 
set of linear equations for the expansion coefficients 

{Gaa — ^a)Ea + ^ GafjEfj — G^E'^ — la, 

(G^^ — Ti^)E'^ + G^^El — G^ E^ = 0. (6) 
The coefficients En and E' are 



■'a I 



Br, 



(7) 



Ea — + Bn 

E'n^- {Aat 

so that the system of Eqs. Q connects the electric 
field modal amplitudes on the incoming interface, z = 
0, and on the outgoing one, z = h. The term G„ 
describes the coupling between the input and output 
sides of the holes, and arises from the reflection of 
the waveguide mode at the openings: 



G 



iYa 



o2iq^ah 
a2iqzah 



(8) 
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where Yq, = qza/g is the admittance for the TE waveg- 
uide mode, and Yq = g/qza is that for the TM one. 
The coupling matrix elements of the system of Eqs. ^ 
are related to the in-plane components of the EM 
Green function dyadic G. The latter is associated to a 
homogeneous medium in three dimensions and is rep- 
resented in the waveguide mode space: 



G^p = {a\G\P) 



'a\K){KY 



(9) 



Here is the admittance of the mode in free space: 
Yk,s = kz/g and lk,p — g/kz- The right-hand side 
(r.h.s.) term takes into account the overlap between 
the incident plane wave and the waveguide mode \a) 
inside the hole. Considering a normalization for the 
incident wave such that the energy flux through the 
hole is unity, Re(/^^jg dS x H*) = 1, we obtain 



/„ = 2iJY, 



(10) 



where Y^; is the admittance of the incident plane wave 
in free space. 

Once the solution of the system of Eqs. ^ is found 
and the modal amplitudes are known, the normalized- 
to-area transmission coefficient can be written as 



T 



E 



lTn{G^p)E'^E'. 



(11) 



Until now we have not mentioned any restrictions 
on the hole shape, which in our formalism only influ- 
ences the structure of the waveguide modes \a). In this 
paper, however, we restrict ourselves to the considera- 
tion of both circular and rectangular holes, where the 
waveguide modes are known analytically. 

For thick films the solution of Eqs. ([6| converges 
quickly. The reason is that in the subwavelength 
limit the amplitudes of the waveguide modes decay in- 
side the hole. As the decay is characterized by the 
propagation constants qa, the higher the waveguide 
mode index, the weaker its influence on the transmis- 
sion. Therefore, only a few waveguide modes (with the 
smallest decrements) contribute the transmission. 

In contrast, when the thickness h of the PEC tends 
to zero, the solution of the system of Eqs. (|6| involves 
many waveguide modes (hundreds or even thousands) 
to provide the precise result. However, we shall show 
below that a very accurate computation of the trans- 
mission can be performed with some tens of waveguide 
modes. We shall find the asymptotic value of the trans- 
mission by a fitting of the convergent result. 

To conclude this section we shall show how to ap- 
ply the above equations when the bounding dielectric 
media have permittivities e/, £///. In this case, for 
a medium with the dielectric constant e, the propaga- 
tion constant of the mode and the admittances become 
kz = i/e^^ _ fc2^ Yi^^s = kz/g and Yk,p = eg/kz- Then 



the wavelength-dependent tensor Gap = G'q^(A) de- 
scribing the interaction of the waveguide modes inside 
the cavity through the EM continuum in vacuum be- 
comes a function of e, Gap{X; e). In the upper equation 
of Q Gai3 changes to GapiX; £/), whereas in the lower 
to Ga§{\ ^iii)- Working out the expression given by 
Eq. (l9|, we find that the relation between the tensors 
is 

Ga0{X;e) = ViGap{X/V~e;e^l). (12) 



III. SMALL HOLE LIMIT 

In this section we show how to simplify computation 
of the transmission when the linear size of the hole, 
a ^ \/S is small compared to A. More precisely, we 
consider the limit 



s = ga ^ 1. 



(13) 



In this limit an accurate numerical computation of the 
tensor Gap becomes problematic. This is related to 
the orders of magnitude difference between imaginary 
and real parts of Gap- We have found that in the low- 
order in parameter e the non- vanishing elements of the 
tensor depend upon e as (see Appendix [A| 



lin{Gap) 



Re{Gap) - l/e. 



(14) 



Both real and imaginary parts of the tensor are impor- 
tant in spite of their substantial difference: Ile{Gap) 
contribute into the amplitude of the waveguide modes 
(see below), while Im(G'Q^) takes into account the ra- 
diation into free-space [see Eq. (111]. For e <C 1 the 
imaginary part of the tensor allows analytical compu- 
tation, and the real part can be considerably simplified 
(see details in Appendix [A| . 

In this paper we restrict ourselves to normal-incident 
wave transmission. In this case the impinging wave 
can only couple to certain waveguide modes of TE 
type. The analysis of the Gap elements shows that 
the contribution to the transmission from TM waveg- 
uide modes is always negligible. This results from both 
a weak coupling between the waveguide modes of dif- 
ferent polarizations inside the small hole, and a weak 
coupling of TM waveguide modes to the EM contin- 
uum of vacuum half-spaces. 

For the circular hole the incident plane wave couples 
directly with only the "horizontal" TEin waveguide 
modes with integer n (where the first index indicates 
the number of semi-periods of the field placed along 
the polar angle). Coupled between themselves, these 
waveguide modes are the only ones contributing into 
the transmission. 

For a rectangular hole, when the electric field is di- 
rected as shown in Fig. [T] the illuminated waveguide 
modes are TEon, with odd n (the first and the second 
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indices define the number of semi-periods of the field 
placed along x and y directions respectively). Only 
this set of waveguide modes appears in the summation 
for the transmission according to Eq. (Ill [all other 
elements of lm{Gai3) are negligible]. However, TEon 
modes couple to TEmn ones with even m and odd n 
through the system of Eqs. ([6]) and must be taken into 
account . 

Due to property expressed in Eq. (14), in the ex- 



treme subwavelength limit the transmission coefficient 
scales as 



T = eV(/i). 



(15) 



The thickness- and shape-dependent function '4}{h) is 
given by the solution of the system of Eqs. ([6]) with 
appropriately normalized coefficients 



where 



V'(/i) = ^G„0^„(/i)^^(/i), 

a, 13 



Gap = ^ 2 StUd Ea — — ~ 



(16) 



(17) 



The amphtudes Ea satisfy the system of Eqs. ([6|, 
where the imaginary part of the Green tensor is ne- 
glected and its real part must be normalized with 
the small parameter Ke{Gap) eKe{Gap)- The 



coefficients of Eq. ^ are replaced by 



sG' 



A. Perfect electric conductor screen 

In this subsection we compare the solution based on 
our formalism with some known results on the trans- 
mission through apertures in an infinitesimally thin 
PEC screen. By reaching an excellent agreement with 
these results, we justify the applicability of the modal 
expansion even in the most unfavorable conditions for 
it. 

In the limit h ^ 0, special care is needed when solv- 
ing system of Eqs. ^ due to the divergency of the 
coefficients and Sq. In order to remove this di- 
vergency, we expand the amplitudes of the waveguide 
modes over gh 



Ea = E^^^ + {gh)E'i^ 4 
E'^^ET^ + {gh)Ei'^ 



(18) 



Equating terms proportional to [gh) ^, we obtain 

i?i°)=-i?f). (19) 



Then, keeping terms of zero-order in gh, we have from 
Eqs. ^ 

Y^GapEf^ ~Ya{E^2~^+E'^^^)=Ia, 

Ga,Ef + YaiEi'^ + E'^^^) = 0. (20) 

/3 



Adding the two equations in (20) and neglecting the 
imaginary part of the tensor Gap, due to Eq. (14), we 
arrive at the final system of equations 



2ymGa0)Ef 







la 



(21) 



Now, when the film has been converted into a screen 
with infinitesimal thickness, the amplitudes of the 
waveguide modes at the incoming and outgoing faces 
of the film are equivalent. Therefore, the terms Eq, 
and G^ responsible for the reflection and coupling of 
the waveguide modes inside the cavity are not present 
in Eq. (21), and the waveguide mode amplitudes are 



coupled by the doubled elements of Re(Ga/3). 

Using the normalization defined by Eqs. (16), (17), 
the transmission coefficient can be written in the form 



(15), where the function Tp{h) becomes a constant de- 



fined by the shape of the hole. 



(22) 



The simplest solution of the system of Eqs. (21 1 is 



obtained by retaining only the fundamental waveguide 
mode. For the circular hole, the Green tensor element 
corresponding to the fundamental waveguide mode is 
GtBhTBh = 1.1951/e 0.2789ie2. For the square 
hole GteoiTEoi = 0.9577/e -I- OMiie^. Within this 
single-mode approximation, taking the r.h.s. from Ap- 



pendix A, the transmission pre-factors of Eq. (22) are 
found immediately: for the circular hole Go = 0.1634 
and for the square one Gq — 0.3041. For a circular 
hole, this minimal model provides the transmission of 
order of 30% with respect to the exact Bethe's result 
Go = 64/(277r2) ~ 0.2402. Therefore, more waveguide 
modes must be taken into account in order to obtain 
the correct result. 

The convergency of the constant for the circle. Go, 
and for the square, Gg, with respect to the number 
of waveguide modes are shown in Fig. [2j After having 
computed the value of this constant for several tens of 
waveguide modes, we fit it by a polynomial 



C 



E 

m=Q 



am 



(23) 



where N is the number of the waveguide modes. 

The value ag gives us the constant G. In the cal- 
culations shown in Fig. [2] the fitting has been done 
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FIG. 2: (Color online) Normalized by (ga)^ transmittance 
for the circular and square apertures of the size a in the 
PEC screen. The asymptotic values are shown by the 
dashed lines. 
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FIG. 3: (Color online) Transmittance for the rectangular 
aperture in the PEC screen, normalized by (ga^)^ (gay)^ . 
The shape of the rectangle for a^/ay equal to 1/3, 1 and 3 
is presented. 



by a forth-order polynomial {m^ax — 4) using 50 
waveguide modes. For the circular hole we obtain 
Co — ao — 0.2403 with an error of only 0.05% with 
respect to the exact Bethe's value. For the square hole 
we obtain Cq ~ 0.4565. Notice that the normalized 
transmittance through a square hole with the side 2a 
is about two times larger than that of the round hole 
with the radius a. 

In the case of rectangular holes, it is useful to write 
the transmittance as T = e^e^C(T), where 



lyQ and T is the aspect ratio 

(Xx 1 0>y . 



(24) 



In this representation C(t) reflects the dependency of 
the transmittance upon the aspect ratio for a constant 
area of the hole. As seen from Fig. [3] this depen- 
dency is a fast function of r. This is due to a strong 
dependency of th e polarizability of the hole upon the 
aspect ratio .'i^'^ From the point of view of the modal 
expansion formalism, the cutoff wavelength of the fun- 
damental waveguide mode is Ac = 4ay, when the field 
is parallel to a;-axis. Therefore, the larger the ay, the 
closer the hole to the resonant regime 

We have fitted the dependency shown in Fig. [3] in 
the region r G [1/3, 3] by the following function 



C(r) = 0.0132 -I- 0.2127/r -I- 0.2174/t2 



(25) 



This fitted function provides an excellent approxima- 
tion to the transmittance: in Fig.[3]the curve given by 
Eq. ( 25 1 is indistinguishable from that obtained from 
the numeric calculations. In the interval r < 1 the 
dependency C(t) can be extracted from Ref. [THl and 
we have checked that it coincides with Eq. (251. 



Let us now turn to the dependency of the transmit- 
tance upon dielectric permittivities of the bounding 
media. Addressing to the general property of Ggp 
in Eq. ( [T2| and to the scaling given by Eq. (14 1, 
where e ~ 1 /A, we see that the imaginary part of 

the tensor is a function of ej, namely Im(G'^^) ef^^, 
where I ~ I, III. In contrast, Re(G'^^) is not depen- 
dent upon £/, so that the amplitudes of the waveg- 
uide modes depend upon e/ only through the r.h.s.. 
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Then it follows directly from 



Eq. (11) that the transmittance for arbitrary substrate 



and superstrate is related to the one when the system 
is in vacuum as 



T(A, ei,eiii) - Jeie]jjT{\ a = 1, em = 1) (26) 



for e^/ei <C 1, s^ftiu <C 1. Note that the transmis- 
sion is not symmetric with respect to the region of in- 
cidence: r(A, 67,6///) 7^ T(A, £///, £/). This may seem 
to be paradoxical, as the transmission coefficient of the 
incident plane wave into the plane wave with the same 
in-plane wavevector is symmetric due to the symme- 
try of the scattering matrix. The total transmittance 
r(A), however, takes into account the transmission of 
a plane wave into a continuum of states. As the den- 
sity of final states depends on the dielectric constant in 
the transmission region, so does r(A). The integration 
over the scattering amplitudes yields a factor propor- 
tional to the wavevector squared modulus in the trans- 
mitted medium, T ^ (^e///g)^ ^ £///, and breaks the 
1-111 symmetry in Eq. ( 26 ) . 
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B. Thick and medium films 

The EM fields inside a subwavelength hole decay ex- 
ponentially with both the film thickness and the prop- 
agation constants of the waveguide modes. In the limit 
of a very thick film, e"^'"^'"''' <C 1, the coefficients of 
Eq. (l8| are simplified and the solution of the system of 
Eqs.T6| can be cast in the matrix form 



E 



D-'I, 



E' 



D-^dD-^L 



where 



D = \\Gafi + i6a,i3Ya\\, d = -2i\\5a,pYae 



(27) 



|.(28) 



Thus, for very thick films the waveguide mode ampli- 
tudes on the input side do not depend upon h, and are 
of the same order as the incident field. Conversely, the 
field amplitudes on the output side are exponentially 
decreased. 

For very thick films only the fundamental TE waveg- 
uide mode is expected to contribute into the cou- 
pling between both sides of the hole. Therefore, we 
can retain in the diagonal matrix d only the element 
corresponding to the fundamental waveguide mode. 
This means that for a thick film the transmittance 
decays as T ~ Q-'^\<izo\h^ where qzo is the propa- 
gating constant of the fundamental waveguide mode. 
For a rectangular hole with the sides 2ax and 2ay, 
Qzo — \/ 9^ ^ ('''/2ai/)^ — «7r/2ay, and for a circular 
hole of the radius a, q^o — \/ 9^ ~ (ui/a)^ ~ iui/a 
(see the definition of Um in Appendix A 1 ) . Then for 



arbitrary film thickness it is useful to rewrite T given 
by Eq. (15 1 in the following form 



y — £4g-2|g^ol'i 



(29) 



The computations of C{h) in the interval h/a G [0, 1] 
for both the square and circular holes are shown in 
Fig. |4] In both cases we phenomenologically adjust 
this dependency by the function 



C{h) = C°° + {C- C°°) e-^''/", 



(30) 



containing fitting parameters C, C°° and 5. In the 
limit of a screen, /i = 0, Eq. (29 1 transforms into 



Eq. (22 1, while for the infinite film thickness f{h) be- 
comes a constant C°° = C{oo). The constants C have 
been given in Subsection |III A[ The values of the con- 
stants for the infinite film thickness have been found 
to be ~ 0.1694 and ~ 0.3027. The values for 
parameters S are So — 6 and (5n ~ 5. 

For the rectangle the constants C and C°° in 
Eq. (30 1 are functions of the aspect ratio. C — C{t) 

^(t) reads 



is given by Eq. ( 25 1 , and C 



C°°{t) = 0.2298/t + 0.08262/r2 



(31) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 08 0.9 1 
film thickness, h/a 

FIG. 4: (Color online) The normalized transmission 
through small circular and square holes as function of the 
film thickness. Calculations using only the fundamental 
waveguide mode are represented by dashed lines. The val- 
ues for the fitted expression (|30[) are marked by " x " . 



Thus, the formula (29| is applicable for the rectangle as 
well, if we replace s 
h/ay in Eq. (30) (we find that 



in Eq. ([29| and h/a 
6 for the rectangular 
hole coincides with that for the square hole 6a). 

In order to make our results more accessible, we sum- 
marize in Table |l] the constants appearing in Eqs. ( p9| 
and ([30|. 

As we see from the dependencies shown in Fig. |4] 
the difference between the full calculation and single- 
mode approximation deceases when h increases. For 
the zero thickness screen it is of order of 30%, and for 
the limit /i ^ oo it is of order of 10%. However, if 
we extract the amplitude of the fundamental waveg- 
uide mode from the full (many-mode) solution, and 
use only this waveguide mode to compute the trans- 
mission according to Eq. (Ill, the difference between 



this calculation and the exact value is reduced. 



TABLE I: Resume of the parameters appearing in the an- 
alytical representation of the transmission, see Eqs. ( |29| , 



Circular hole Square hole Rectangular hole 
(radius=a) (side=2a) (sides=2a2,, 2ay, r = ^ 



C 



C° 



(Bethe^) 0.4565 0.0132 



0.2127 I 0.2174 



0.1694 



0.3027 



S 6 



2a„ ' 



To summarize this section, in the extreme subwave- 
length regime, many waveguide modes are necessary to 
provide the precise value for the transmittance. On the 
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other hand, the fundamental waveguide mode plays a 
crucial role in the process, especially in the coupling 
between the fields on top and bottom faces of the hole. 
As we will show, this reflects the fact that the funda- 
mental waveguide mode possesses the largest induced 
dipolc moment (see Section |v] below) . 



IV. HOLES OF A MODERATE SIZE 

The sizes of apertures used in the experimental 
samples are often not in the extreme subwavelength 
limit considered in the previous section. For ex- 
ample, in original experiments on enhanced optical 
transmissiorPl the sizes of the holes were of order of 
200 — 300 nm, so the parameter e was e > 1 in the 
visible range of the spectra. 

In this section we study holes of moderate sizes: still 
in the subwavelength limit, but with the condition ( 13 1 
not fulfilled. However, we still consider wavelengths 
larger than the resonant wavelength of the hole. For 



the resonant transmission through a single hole at 
wavelengthes close to the cutoff, we refer the reader 
to Ref. TH1: 

When the size of the hole increases, the transmit- 
tance value can be refined by retaining the next terms 
in the expansion over e. For a circular aperture in the 
PEC screen of zero thickness a few first terms were 
computed by Bouwkamp.'^ However, such a series has 
radius of convergence i?e ~ 1, and therefore is appli- 
cable in a narrow region of e. 

For arbitrary film thickness and size of the hole, the 
transmission can be accurately computed numerically 
using the modal expansion. But from the practical 
point of view it is useful to have an analytical formula 
containing the parameters of the hole. We have nu- 
merically computed the transmission through a square 
hole retaining many waveguide modes in the system of 
Eqs. ([g]) in a wide range of the size 2a, and the thick- 
ness h. Then we have approximated the calculations 
generalizing dependency ( 29 1 by adding a quadratic in 
£ term to the function C(Ti) 



T - 



(32) 



where C2{h) has the same form that C{h) in Eq. (30 1 
C2{h) = + (C2 - C^) e'^^l". (33) 

The adjusting constants are C-iu — 0-66 and ~ 
0.43. This dependency provides the transmission with 
an error not exceeding a few percents in the region 
a/A < 1/6 (i.e. e < 7r/3), and for an arbitrary film 
thickness h. In the limit of the zero thickness screen, 
Eq. (32 1 has a form similar to Bouwkamp's expansion 
(up to the sixth-order term) However, C2 is not a 
constant defining the sixth-order term of the authentic 




FIG. 5: (Color online) The normalized transmission 
through small rectangular hole, C(K) +£^Cq{h) as a func- 
tion of the film thickness, h/a and the wavelength. A/a, in 
the units of the hole half-side. 



expansion over e. The latter term in the real expansion 
only refines the value of the transmittance in the small 
hole limit, meanwhile C2 results from the fitting of the 
spectrum captured in a wider region. 



V. THE HOLE AS A MULTIPOLE 

It is well known that the EM field of any source 
can be considered as resulting from a superposition 
of multipoles (see e.g. Ref. H]and references therein). 
Recently this viewpoint has been used to study the 
transmission properties of collections of holes.'^Sl In 
this section we establish the connection between the 
mode-matching formalism and the multipole expan- 
sion. For this purpose let us focus our attention on 
the far-field. The in-plane components of the electric 
field are given by Eqs. (|3]), and its z-component can 
be derived from Maxwell's equations by a straightfor- 
ward differentiation. The expression for the far-field 
can be derived with the help of the scalar free-space 
Green function associated to the Helmholtz equation 
in three-dimensions.'^ Additionally, in the region of 
transmission, the expression for the E far field, can be 
computed bu using Green's function identities.^ The 
expression for E in terms of the integration of the field 
at the face of the hole Et(ri) over the hole area is 



27r 



u X dS'nx E;(r;)e-'»"""', (34) 



where u is the unit vector pointing into the observation 
point, u = r/r, and n is the external normal. Eq. (34l 
has the form of retarding potential^ resulting from the 
charges induced by the electric field on the face of the 
aperture. The field an the face of the hole z = h, is 
expressed through the waveguide mode amplitudes 



(35) 



In the reflection region, I, the scattered far-field has 
a form similar to Eq. (34 1, but in terms of the field 
on the interface z = 0, i.e., —E'^ is replaed by in 



Eq. (351 



If we expand the exponent in the integral of Eq. ( 34 1 , 
g-igur _ 1— igur'+..., the far-field can then be written 
in the form of effective multipole^ 



E^"''(r) = g^ — (u X p -Km - 



u • Qm + ...) X u, 



(36) 



where p and m are the effective electric and magnetic 
dipole moments respectively, Qm is the effective mag- 
netic quadrupole tensor, etc. Comparing Eq. (36) with 



Eqs. (34) and (35 1, we conclude that each waveguide 



mode can be interpreted as a superposition of effective 
multipoles; the far field then results from the contri- 
bution of the multipoles of all the waveguide modes. 
For example, the effective magnetic and electric dipole 
moments of the waveguide mode \a) read 

= J dS' nx (r'tla), 
s 

Po.^^JdS'r[x[nx{r[\a)]. (37) 

s 

Denoting the dipole moments of the reflection region 
by p and m, and those of the transmission region by 
p' and m', we write them as 

a a 

' = -5^i?>„, p' = -^i?;,p„. (38) 



m 



m 



Thus the hole of a flnite thickness can be considered as 
a coupler between induced multipole moments at both 
faces. For example the systems of Eqs. (38) and ^ 
describe the coupling between dipole moments that is 
enough in the low order approximation, when the hole 
is small. We would like to stress that the effective 
dipole moments are coupled through all the waveguide 
modes presented inside the hole. For thick films, how- 
ever, the contribution of the fundamental waveguide 
mode dominates over the other modes. 
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(a) (b) (c) 
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FIG. 6: (Color online) The magnetic dipole moments and 
the electric field amplitude spatial distribution (at z = h/2) 
for the square and circular holes. The amplitudes of the 
dipole moments are shown by the bars. The radius of the 
circular hole is a, and the side of the square hole is 2a. The 
considered PEC slab has thickness h = 0.1a. The distri- 
bution of the total electric field modulus together with the 
electric field lines are shown in (a) and (d). The contribu- 
tion into magnetic moments from the two lowest waveguide 
modes and the fields of these waveguide modes are also 
shown for both hole shapes: TEoi and TEos of the square 
are represented in (b) and (c); TEn and TEis of the circle 
are shown in (e) and (f). 



Only waveguide modes with zero effective electric 
dipole moment couple to normal-incident light and 
generate the far-field. For the circular hole calcula- 
tions using Eq. ( 37 1 show that only "horizontal" T£'i„ 



waveguide modes with integer n contribute to the mag- 
netic dipole moment. For the rectangular hole the 
contributing waveguide modes are TEq^ and TE^q 
with odd m. When the incident electric field is di- 
rected along Ox (magnetic field along Oy), the mag- 
netic dipole moment induced on the incoming face of 
the hole with the area S i^^U 



m — fiSey, where M = / E^Ca 



(39) 



We obtain that for the circular and rectangular holes, 
the constants C„ have the values 



I 



*v/2vr3(u2-l)^ 



Greet 



2t 



(40) 
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where n is integer and m is odd. For the dipole mo- 
ment m' induced on the outgoing face of the film 
must be changed to —E'^. We have checked that the 
value for ^ of the circular hole in the limit of zero- 
thickness screen (that we obtain with the help of 50 
modes), — 0.07624, reproduces the Bethe's value 
4/(37r5/2) ~ 0.07622. For the square hole we have 
found /in = 0.08253. 

Fig. [6] renders the effective dipole moments and the 
exact fields of the small circular and rectangular hole 
in the PEC slab of the finite width. The contribution 
from the fundamental waveguide mode is dominant, 
so that the field distribution inside the hole is very 
similar to that of the fundamental waveguide mode. 
Interestingly, the vector lines of the full field are similar 
for square and circular holes. 

It must be noted that in our method the effective 
dipoles must be computed after solving the systems of 
Eqs. (|6| , ( 38 1 . We have been unable to derive the equa- 
tions governing the effective dipoles directly, without 
previous calculation of the amplitudes Ea, E'^. 



VI. CONCLUSIONS 

To conclude, this paper has explored the EM trans- 
mission through both small and medium-size isolated 
holes in a perfect electric conductor screen of arbi- 
trary thickness. We have used the modal expansion 
and have shown that this technique is applicable even 
in the limit of zero film thickness. This latter limit 
has been used to check the correctness of the theory 
through the comparison with Bethe's result. We have 
phenomenologically fitted the transmittance in a wide 
region of parameters by simple analytical functions. 

We have connected the formalisms based on either 
modal or multipole expansions, and showed that the 
induced dipole moments are coupled via all the waveg- 
uide modes inside the hole. Our results indicate that 
the fundamental (lowest-order) waveguide mode pos- 
sesses the largest dipole moment. It is also responsible 
for the attenuation of the transmission with the in- 
crease of the film thickness, T ~ e"^'"^^"'''. 
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APPENDIX A: COMPUTATION OF THE 
GREEN'S TENSOR 

Here we give explicit expressions for the tensor Gap 
in the small hole limit. When the parameter e is small, 
the integral in Eq. (|9| can considerably be simplified. 
We illustrate these simplifications on the example of 
the diagonal element for the fundamental mode of the 
square hole with side 2a. This diagonal element has 
the following form in polar coordinates (q^ — qcos9, 
qy gsinfl, q = k/5) 



t/2 



GtEoiTEoi = 8i / d9 dq 



1 - sin^ ( 



q^T^ 
00 

[1 — cos(2egcos0)][l -I- cos(2e(7sin 0)] 



(Al) 



We have taken into account the parity of the integrant 
both in qx and qy, and, therefore, integrated over the 
first quadrant of the q^-qy plane only. 

We see that the integrant in Eq. ( Al I is either purely 



real or purely imaginary depending only upon the 
square root \f\--^ . Therefore, the integral over q can 
be separated into two integrals, one from to 1 (yield- 
ing the imaginary part of GtEoiTEqi) and the other 
from 1 to CX3 (yielding the real part of GteqiTEoi)- 
Then we expand the integrant for the imaginary part 
into a series over the parameter e, retain the leading- 
terms only and take the integral analytically. The re- 
sult is lm{GTEoiTEoi) - 32eV(37r3). To treat the in- 
tegral for Re(G'T_EoiT£;oi)' make the change of the 
variable ^ = 2q£. Taking into account that the region 
^ < 2e weakly contribute to the integral, we simplify 
the square root as \/\ — q^ ~ iq. Then we replace 
the lower limit ^ = 2e for ^ = and approximate 
1 — Qje^ by —^sva^Qje^ in the denominator, 

so that the real part of GTE^TEm becomes 




[1 - cos(^ COS 9)] [1 -I- cos(^ sin 0)] 

(^2sin2 6'-7r2)^ 



X 



(A2) 



In order to perform the integration over f analytically, 
we extend the integrant into the complex plane, chang- 
ing the trigonometric functions of 1^ to the exponential 
functions. Finally, Eq. (A2| is derived applying the 



Residue theorem, by taking into account the presence 
of the poles on the real axis ^. The integral over is 
taken numerically. 

The other non-vanishing tensor elements both for 
circular and rectangular holes are simplified analo- 
gously. Note, however, that for a circular hole it is 
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more convenient to perform the analytical integration 
over 9 in the real part of Gap, using the identities for 
the Bessel functions. The integral over ^ can then be 
performed numerically. 

In the two following subsections we give simplified 
Green's tensor elements for the rectangular and circu- 
lar holes. 



The right-hand side term of Eq. ^ is 



(A5) 



1. Circular hole 

Consider a circular hole of radius a. The imaginary 
part of the Green tensor is approximated as 



Iva^GTEir^TE^^,) — 



(A3) 



where are the solutions of the equation J{(w„) = 0. 
After performing the integration over 9, the real part 
reads 

o 

Re(GT_Bi„T£;i„, ) 



2. Rectangular hole 

For a rectangular hole with the sides and 2ay 
the imaginary part of Gap simplifies to 



Im(G' 



TE„^^TE„,,) = Sn oSn'.OT, i (^6) 



where ex,y = o,x,yg- The real part takes the following 
form 



Tr/2 oo 



(A4) where 



nm:n m' 



2CT„cr„' y/{n^ + m^T'^){n'^ + m'^r^)^^ [1 - cos(C cos 9)] [I + cos(^r-i sin 9)] 



sin^ 9 cos^ 9 



\cose) ? \cose) ? V sin 9 / ? \ sin 9 J 



(A8) 



If rt = then (T„ = \/2 a nd (7 „ = 2 otherwise. To sim- Eq. ([6| 
plify the double integral ( A7 ) we transform the cosines 



IS 



in the nominator of <i> into exponential functions. The 



integration over ^ is performed analytically, extending Ite„^ = Sn,o • (^^) 

the integrand into the upper complex half- plane, and ^''^ 
using the Residue theorem. The illumination term of 
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